Chemical sensing by cell-surface chemoreceptor arrays: the roles of receptor 

cooperativity and adaptation 
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Most sensory cells use cross-membrane chemoreceptors to detect chemical signals in the environ- 
ment. The biochemical properties and spatial organization of chemoreceptors play important roles 
in achieving and maintaining sensitivity and accuracy of chemical sensing. Here we investigate the 
effects of receptor cooperativity and adaptation on the limits of gradient sensing. We study a single 
cell with aggregated chemoreceptor arrays on the cell surface and derive general formula to the 
limits for gradient sensing from the uncertainty of instantaneous receptor activity. In comparison 
to independent receptors, we find that cooperativity by non-adaptative receptors could significantly 
lower the sensing limit in a chemical concentration range determined by the biochemical proper- 
ties of ligand-receptor binding and ligand-induced receptor activity. Cooperativity by adaptative 
receptors are beneficial to gradient sensing within a broad range of background concentrations. Our 
results also show that isotropic receptor aggregate layout on the cell surface represents an optimal 
configuration to gradient sensing. 

PACS numbers: 87.16.dr, 87.17.Jj, 87.18.Tt 
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I. INTRODUCTION 

Cellular sensory systems can detect temporal and spa- 
tial changes of environmental signals. For example, bac- 
teria sense chemical gradient by motion that translates 
spatial chemical concentration asymmetry into tempo- 
ral asymmetry. Larger chemotactic eukaryotes are able 
to sense the spatial gradients of chemoattractants across 
the cell dimension. A primary task for a sensory cell is 
to respond to small chemical concentration changes or 
shallow chemical gradients with sufficient accuracy un- 
der stochastic noises. Berg and Purcell, in their classic 
study Pi, showed that the fundamental physical limit of 
concentration sensing is set by the uncertainty of ligand 
diffusion to a sensory cell regardless of biochemical details 
in a sensing mechanism employed by the cell. This result 
was generalized by Bialek and Setayeshgar Q using the 
fluctuation-dissipation theorem, and was also applicable 
to gradient sensing Q. 

Ligand detection by membrane-bound chemoreceptors 
is the first step for chemical sensing in a vast majority 
of sensory systems of living cells. Uncertainty of ligand- 
receptor binding and stochastic dynamics of downstream 
cellular signaling may introduce additional noises. Re- 
ducing such noise may help a system to operate near 
the fundamental Berg-Purcell limit. However, the exact 
sensing limit under proper ligand-receptor interaction re- 
mains elusive. In particular, the effect of receptor cooper- 
ativity is controversial 043] • Receptor cooperativity can 
sensitize the receptor response to small signal changes, 
but in the meanwhile it also amplifies stochastic fluc- 
tuations in signal. Using the Monod-Wyman-Changeux 
(MWC) model, Bialek and Setayeshgar [J] showed that 
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receptor cooperativity helps to lower the threshold of 
concentration sensing to approach the Berg-Purcell limit 
set by ligand diffusion. Hu et ah by an Ising-type 
model, showed that receptor cooperativity improves gra- 
dient sensing within a shortened dynamic range of back- 
ground concentrations. Most recently, Skoge et ah Q 
studied chemoreceptor activity by a dynamic Ising model 
and showed that receptor cooperativity could slow down 
receptor activity and thus may not improve the signal- 
to-noise ratio due to reduced response time. 

Another crucial mechanism found in many cellular sen- 
sory systems is adaptation, which maintains the sen- 
sitivity of a system to varying levels of environmental 
signals. It is known that adaptation tunes kinetics of 
chemoreceptors at the molecular level. For the exam- 
ple of the better understood E. coli, the activity of re- 
ceptors are modulated by multisite receptor methylation 
and demethylation. However, the effect of receptor coop- 
erativity to chemical sensing limit under the context of 
receptor adaptation is yet to be examined. 

Using a theoretical model, here we study the role of 
receptor cooperativity on the accuracy of chemical sens- 
ing under the influence of receptor adaptation. We con- 
sider a sensory cell with a hierarchical organization of 
cell-surface chemoreceptors. Receptor cooperativity is 
described by the classic MWC model [1] , which was orig- 
inally developed to explain allosteric regulation of multi- 
subunit proteins and has been widely used to model 
chemoreceptor coupling in bacteria [9Hl3|. We derive 
formula for the limit of gradient sensing based on the 
uncertainty of instantaneous receptor states at the equi- 
librium. Our results show that cooperativity by non- 
adaptative receptors reduces instantaneous noise of re- 
ceptor states within a limited range of background con- 
centrations determined by the biochemical parameters 
of receptor dynamics and ligand-receptor binding. In 
contrast, cooperativity by adaptative receptors improves 
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the sensing accuracy across a wide dynamic range. We 
also show that the layout of receptor aggregates on the 
cell surface significantly affects the sensing limit with the 
isotropic layout being the optimal. Although the sensing 
limit is sensitive to the cell orientation under anisotropic 
aggregate layouts, the effect of receptor cooperativity is 
invariant. 



II. THEORY 
A. A two-dimensional cell model 

We consider a model (Fig. [IJ for a two-dimensional 
chemotactic cell subject to a chemoattractant gradient 
field C?(x;p, </>), where x is the spatial coordinate, and 
p and 4> are the steepness and direction of the gradi- 
ent. The steepness is defined in a polar coordinate sys- 
tem as: p = —4 s , which is a normalized concentration 
change across a reference distance ro along the direction 
of the gradient. The cell has M receptor aggregates dis- 
tributed at distinct locations on the cell surface. Coop- 
erative receptors form clusters of size n and each aggre- 
gate contains a number of such independent and non- 
interacting receptor clusters. Assuming aggregate to has 
N m receptor clusters, we have the total number of re- 
ceptor monomers in a cell as: A to t = n ^2m=i N m - This 
hierarchical organization of receptors on the cell mem- 
brane is general and can be parameterized to study a cell 
that has aggregated receptor arrays with (n > 1) or with- 
out cooperativity (n = 1) or has non-aggregated receptor 
monomers (N m — 1 and n = 1). 

In the absence of ligand, a receptor switches between 
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active ( "on" ) and inactive ( "off" ) state with a free energy 
difference AE = E Q {f — E on (in units of kgT, where kg 
is the Boltzmann constant and T is the absolute tem- 
perature). As in an MWC model, coupled receptors in a 
cluster switch in an all-or-none fashion between the "on" 
and "off" states. A ligand binds to a receptor of the "on" 
or "off" state with a dissociation constant K on or K a s, 
respectively. Ligand binding will shift the free energy dif- 
ference between the two receptor states and such change 
can be sensed by a cell to measure the ligand concen- 
tration. At the equilibrium, a receptor cluster to at the 
"on" or "off" state bound to r ligands has the free energy 

of nE on - In (^Jf-) or nE oS - In , r = 0, n, 

respectively, with a multiplicity of (™) . From Boltzmann 
distribution, one can obtain the equilibrium probability 
for a receptor cluster in aggregate to to be active is given 
by ' 
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where c m = [L] m /K s is a normalized ligand concen- 
tration. The ligand concentration [L] m at the location of 
aggregate m is determined by the gradient field Q(x; p, c/>). 
The non-dimensional coefficient a = K a ff/K on . 

At the equilibrium, an instantaneous configuration of 
receptor cluster states measures the chemical concentra- 
tions around the cell. Such configuration fluctuates in 
time due to randomness of ligand-receptor binding and 
stochasticity in receptor state switching, which underlies 
the uncertainty in gradient sensing. Based on the above 
model, below we shall derive the best achievable limit for 
the sensing uncertainty. 

Under the chemical gradient field Q(x;p, <ft), the log 
likelihood for observing a specific configuration of cell- 
surface receptor states is: 
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where k m is the number of active receptor clusters in ag- 
gregate to. An efficient unbiased estimator to parameters 
associated with the gradient field Q (x; p, <fi) has a variance 
limited by the optimal Cramer-Rao lower bound (CRLB) 
that can be determined from the likelihood function of 
Eq.@. For convenience, define (^1,^2) =Jp, We can 
calculate the Fisher information matrix 11 



FIG. 1. (color online). A 2-D chemotactic cell with M 
chemoreceptor aggregates (patches labeled with polar coordi- 
nates (r m , 9 m )) distributed on the cell surface (enclosed curve) 
under a gradient field (J(x;p, <^>) indicated by the arrow. Each 
receptor aggregate has N m ,m — 1,...,M, cooperative recep- 
tor clusters of identical size n. 
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where the expectation (■) is taken with respect to the 
binomial distribution f(k m ;P m: N m ) = (^ m )P, k n m (l - 
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P m )N m -k m _ We then have 
/ M M , 



One can now obtain the Fisher information matrix: 



N m Pm dP m h - N t P[ dPl 



i iri p ™( 1 ~ p ™) ~ P; ) I 

= ^ ^ (fan - N m P m )(k t - NtPj)) dP m dPj 

fki hi Pm ^ ~ (1 - Pi) d^i d^j ' 

Notice that the covariance ((k m — N m P m )(ki — NiPi)) = 
when to ^ I due to independence of receptor aggregates. 
For to — I, the variance of k m is ((fc m — N m P m ) 2 ) = 
N m P m (l — P m ). Therefore, the Fisher information ma- 
trix is: 
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where one can verify that 

dP m „ „ * I - a 
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We define the coefficient 
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which is a function of local chemical concentration c m , 
size of the receptor aggregate N m , the strength of recep- 
tor coupling n and receptor affinities to the ligand K on 
and K tf . We first notice that for the single receptor ag- 
gregate m sensing its local chemical concentration c m , 
the Cramer-Rao lower bound of the sensing variance is 
given as: 



J_ = (l + c m ) 2 (l + ac m ) 2 
uj m N m n 2 P m {\- P m )(l-af 
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The noise-to-signal ratio is determined by a 2 m /c m . 

For chemical gradient detection by all cell-surface re- 
ceptor aggregates, the sensing limit to steepness p or di- 
rection <f> is set by the CRLB's (diagonal entries of the 
inverted Fisher information matrix): 
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To obtain analytical results, we assume that the cell 
resides in a linear gradient field. Across a typical cell 
size (e.g., about 10 /im for Dictyostelium and 1 fim for E. 
coli), the linear gradient is a reasonable approximation. 
c m at the polar coordinate (r m ,6 m ) is: 



co[l + f3 m p cos(0 ■n 
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where cq is the background concentration at the origin 
and the coefficient /3 m = r m /ro is the distance from ag- 
gregate to to the origin normalized by the reference dis- 
tance r . We have 
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In all cases, a 2 and cr 2 , can be numerically computed, 
and it is possible to obtain their analytical solutions un- 
der special cell geometry and surface layout of receptor 
aggregates. Here we specialize to a circular cell of ra- 
dius r and designate the origin as the cell center such 
that j3 m — (3 — r/ro, for all to. For a shallow gradient 
across the cell length (p/3 <C 1), we approximate ui m « ujq 
for all receptor aggregates. lu is calculated by Eq. © 
evaluated at concentration Co, which can be considered 
as the average ligand concentration at the cell location 
(assigned to the cell center). For an equidistant layout 
of identically-sized receptor aggregates over the cell sur- 
face (an isotropic distribution), we obtain the CRLB's 
for uncertainties in sensing p and <f> [l4j : 
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The factor £„ = nPo(l — Po) is related to receptor coop- 
er ativity n, where Po is calculated by Eq. ([T]) at cq. The 
above limits in Eq.([T2]) have a few properties: (i) Both 
a 2 and are independent of the number of receptor ag- 
gregates M and are insensitive to the gradient direction 
(j>. (ii) Directional sensing improves (cr 2 , decreases) as 
the steepness p increases. Fortuitously, <r? is the noise- 
to-signal ratio of steepness sensing, (hi) The receptor 
system loses sensing capability when ligand binding does 
not differentiate the two receptor states (a « 1). (iv) a 2 
is very sensitive to the cell size via the term N tot j3, where 
iVtot may vary with the cell size by a certain relationship 
and therefore a 2 could scale more strongly than the in- 
verse of normalized cell radius f3 2 . Properties (iii) and 
(iv) hold even though the layout of receptor aggregates 
on the cell surface is anisotropic. 



B. Extension to a three-dimensional cell model 

The 2-dimensional model can be extended to a 3- 
dimensional cell to approach a more geometrically re- 
alistic scenario. In the spherical coordinate system, a 
point has its spatial coordinate (p,"d,cp), where p is the 
radial length, $ and ip are the polar angle and the az- 
imuthal angle. Assume a linear chemical gradient with 
a steepness p and a direction along the unit vector 
g = [cos fig , sin 'dg cos tp g , sin ?? g sin ip g ] . For a receptor ag- 
gregate to at cell-surface location A m , the vector along 
the direction from origin O at the cell center to A m , 
OA m , is a m = r m [cosi? m ,sini? m cosyi m ,sini? m sin^ m ]. 
Concentration c m can be calculated as: 

Cm = co(H a m ■ g) = co(l + /3 m pcos7 m ) , (13) 
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FIG. 2. (color online), (a). Function (1 - 2P )ln(l/P - 1). 
Inequality Eq. (O holds when P £ (P_,P+). (b) Solid 
curves show AP as a function of Co at a = 0.1, AP = 
ln( 1 1 + c e " ) - ln(l/P+ - l)/n (lower two curves) and AP = 



■ l + aco ' 



ln(l/P_ — (upper two curves), for two co- 

2 (□) and n = 6 (o). Upper and lower 



l + CKCQ ' 

operativity levels: n 
curves of same n enclose the regime that receptor cooperativ- 
ity (up to n) helps to reduce the variance Up. The beneficial 
ranges of Co for three different AP values (2.0, 1.0, and 0) are 
marked (by dotted lines spanning boundaries specified by Eq. 
(|20|) '). The dashed line indicates the AP adjusted by precise 
adaptation, in which the adapted level of receptor activity is 
a constant P a = 1/2. 



where j m is the angle between a m and g. The dot prod- 
uct a m • g = r m cos7 m projects a m onto the gradient 
direction g, with cos7 m = sin-i9 m sin^ g cos(ip m — ip g ) + 
cos dm cos fig . We have: 
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c Q /3. m cos j m (14) 
c p m p(siwd m cos"d g cos(</? m - tp g ) 

— COS$ m SUl$ g ) (15) 

c /3 m p sin -d m sin fi g sin(y> TO - ip g ) . (16) 



The Fisher information matrix I can be constructed by 
Eq-Q, with the coefficient uj m obtained by Eq.©. Ana- 
lytical formula like Eq. (Il~2")) can be obtained under special 
cell geometry and receptor aggregate layout on the cell 
surface. We will present our results for a 2D cell and 
expect that general conclusions are applicable to a 3D 
cell. 



III. RESULTS 

Here we investigate the role of receptor cooperativity 
on the gradient sensing limit described by Eq. (|12p . Under 
a fixed number of receptors (N tot ), the effect of cooper- 
ativity can be detected as the direction of change in 
with regard to the receptor cluster size n: 
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Beneficial effect of receptor cooperativity to the sensing 
limit requires d^n/dn > 0, which together with Pq = 
[1 + (e^^yr^-)™] -1 implies the constraint: 
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(l-2P )ln( — -1) <1, 



or equivalently, 
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where P_ « 0.176 and P+ = 1 - P- ss 0.824 are the 
two solutions to: (1 - 2P ) ln(l/P - 1) = 1 (see FigC^a) 
for illustration). The feasible range of AE, a and the 
background concentration c is: 



ln(l/P+-l) 
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within which receptor cooperativity up to the magni- 
tude of n can improve the accuracy of chemical sensing. 
Fig. [21b) shows that the beneficial region of free energy 
difference AE is confined in a banded area from around 
at low ligand concentration to around ln(l/a) at high 
ligand concentration. The bandwidth of AE is a con- 
stant - ln(P+/P_), limited by the strength of receptor 
cooperativity n. A stronger cooperativity results in a 
narrower band. For each given AE, the beneficial range 
of background concentration cq is determined by Eq. (|20[) 
(see Fig^b) for illustration). Results in Figl^b) were 
obtained at a < 1, in which ligand binding favors the 
"off" state of the receptor. The results at a > 1 (not 
shown) are clearly symmetric to those at a < 1. We 
note that the above derivation is under the assumption 
that the receptor system has a fixed total number of re- 
ceptors Ntot- The general conclusion also applies to a 
system that has fixed numbers of receptor clusters (N m ) 
in receptor aggregates [l5j . 

In the following, we analyze the influence of recep- 
tor cooperativity on the chemical sensing limit under the 
context of receptors with or without adaptation. 



A. Non-adaptative receptors 

Without receptor adaptation, the equilibrium level of 
receptor activity changes with the ligand concentration 
Co according to Eq. ((TJ) , where free energy difference AE 
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FIG. 3. (color online). <r| under varying background con- 
centrations. Receptor aggregates have an equidistant lay- 
out on the cell surface, (a)-(d) were plotted with different 
AE (values indicated in plots) under different cluster size: 
n — 1 (non-cooperative), 2 (receptor dimer), 'o'; 6, '*' 
(corresponding to "trimer-of-dimers" in bacteria); and 10, 
(corresponding to an estimate by Ref. (Tlj). o"p were calcu- 
lated by Eq. (|8| (markers) and by Eq. ([12} (solid lines) with 
JVtot = 80, 000, M = 3,p = 0.1,0 = 0, a = O.f and = 1. a% 
behaves similarly and is not shown. 



and binding affinities K on and K a g are unmodulated by 
the receptor activity. Fig. [3] shows that in all cases a 2 as- 
sumes a valley-shaped relationship with the background 
concentration cq. er 2 attains a minimum at a cq that 
satisfies the condition da 2 /dco = 0, i.e.: 



2ac 2 l + n(l - a)(l - 2P )c Q -2 = 



(21) 



At the regime AE < 0, where an unliganded receptor 
is biased to the "off" state, increasing receptor cooper- 
ativity always reduces sensing accuracy [Fig. 0(a)] when 
ligand binding as well favors the "off" state (i.e., a < 1). 
According to Fig. [2jb), AE = — 1 is out of the beneficial 
band at any background concentration cq. Intuitively, in 
this case, receptor arrays have limited capacity (most re- 
ceptors are already "off" even before the introduction of 
ligands) to respond to ligand binding by further switch- 
ing off receptor activity. 

At the regime AE > 0, receptor cooperativity reduces 
a 2 within a concentration range around the optimal Cq 
[Fig. Bfb)]. Such advantage is confined in a narrowed 
dynamic range around intermediate background concen- 
trations [the stronger the cooperativity, the shorter the 
range of improvement] instead of being effective at lower 
concentrations, a desirable region for chemotactic re- 
sponses. The results shown in Fig. G^b) approximate 
those by the Ising-type model of Hu et al. 5] in which a 
receptor maintains the "on" state in the absence of lig- 
and and is only switched off by ligand binding. In fact, 
the model by Hu et al. [5| can be considered as a limiting 



case to our model, by requiring e~ AE <C 1 and a<l. 

By contrast, receptor cooperativity improves sensing 
accuracy below the optimal Co near AE w [Fig. [He)]. 
Fig. EJd) shows that at low background concentration 
(co = 0.1) the system achieves optimal sensing accuracy 
at free energy difference AE w of any cooperativity 
level and receptor cooperativity further decreases a 2 (an 
order of magnitude improvement at n = 10 over indepen- 
dent receptors). We can relate this result to E. coli, in 
which Co = 0.1 corresponds to a background concentra- 
tion of 2 nM MeAsp, a threshold level, at K oS = 20 nM 
and cooperativity predicted about n = 10 as in Ref. [111 ]. 



B. Adaptative receptors 

Near precise adaptation was found in bacteria such 
as E. coli 16]. A ligand concentration change triggers 
a transient response in receptor activity followed by a 
slow decay back to the steady state about the prestim- 
ulus level. In other words, receptor adaptation desen- 
sitizes the steady-state activity to ligand concentration, 
allowing the cell to be responsive to environmental sig- 
nals within a wide dynamic range. Here we show that 
adaptation also conditions the receptor system to allow 
cooperativity to improve signal-to-noise ratio in chemical 
sensing under a wide range of background concentrations. 

We assume that the receptor system is adapted to its 
chemical environment before an onset of change in ligand 
concentration. We also assume that the adaptation ma- 
chinery adjusts the free energy difference AE between the 
two states of the receptor and keeps binding affinities K on 
and K g invariant. For the example of E. coli, receptor 
activity-controlled receptor methylation and demethyla- 
tion act as a feedback mechanism that can adjust the free 
energy gap AE. For precise adaptation, equilibrium re- 
ceptor activity is tuned to be a constant, P a , independent 
of the background concentration. By Eq. ([1]), at aggre- 
gate m, we have the adjusted free energy difference: 



AE = In 



1 



n \P n 



(22) 



Concentration changes around c m induces transient re- 
ceptor activity P m away from P a , which is described by 
Eq.([T]) with the above adjusted AE. 

The optimal sensing achieves at the half-activation 
level P a = 1/2 that maximizes £ n , where = 8(1 + 
c ) 2 (l + ac ) 2 / (N tot nP 2 (1 — a) 2 cl). a 2 attains a min- 
imum of al min = 8(a 1 / 2 + l) 4 /(A^totn/? 2 (l - a) 2 ) at 
the background concentration cq = a -1 / 2 according to 
Eq. ([2l"|) [i.e., [L]q = (X n^off) 1//2 is the geometric mean 
of dissociation constants]. Fig[2Jb) shows (dashed line) 
the adjusted AE to adapted receptor activity P a = 1/2 
as a function of ligand concentration. The adaptation 
tunes AE into the middle of the beneficial region band 
and allows the receptor cooperativity to improve sensing 
of small concentration changes. 
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The adaptation mechanism may also likely maintain a 
constant steady-state activity by adjusting ligand bind- 
ing affinities K an and K g, or adjusting all three ther- 
modynamics parameters altogether (l7j . How adaptation 
tunes these parameters altogether remains unclear even 
for the well-studied E. coli chemoreceptors. Nonethe- 
less, as our model shows, as long as a cell achieves 
an adaptation level (precise or imprecise) that satisfies 
P_ < P a < P+, receptor cooperativity is advantageous 
for sensing small changes in chemical concentration or 
gradient. It has been suggested a precise adaptation of 
P a « 1/3 in E. coli [l8[ in the range specified by Eq. ([19]). 
supporting the idea that receptor adaptation and coop- 
erativity act in concert to improve chemical sensing. 



Effects of anisotropy in the layout of receptor 
aggregates 



The analytical formula in Eq. (|T2j) do not apply to the 
scenario where receptor aggregates anisotropically locate 
on the cell surface. As we will show, such an anisotropy 
may cause substantial dependence of and on the 
gradient direction <f>, the number and locations of recep- 
tor aggregates. We define a metric, TZ e , to quantify the 
layout anisotropy: 



Tie 



2tt/M| 



4tt(M-1)/M 



which is a normalized mean angular dispersion between 
two immediately adjacent aggregates from the average 
angular distance, 2ir/M. Without loss of generality, the 
angles of receptor aggregates are labeled in an order such 
that i+ i > 6i, for i — 1, .., M, with the periodic bound- 
ary condition: 6m+i = 0\ +2tt (see FigfT]). lZ e ranges 
from for the isotropic case to 1 when all receptor ag- 
gregates gather at a single location. 

Each effective receptor aggregate acts as a local con- 
centration sensor, and a pair of aggregates can detect the 
projection of the gradient onto the direction connecting 
the two aggregates. Two independent projections are 
needed in minimum to reconstruct the gradient. In our 
model, the minimal number of effective cell-surface re- 
ceptor aggregates is M = 2 because of the implicit as- 
sumption that one parameter of the linear gradient, the 
background concentration Co, is known. This assumption 
is justified by considering that the cell has a memory for 
a recent history of its chemical environment. 

A sensing singularity occurs when the cell (1) has only 
one effective receptor aggregate because two aggregates 
locate too close to each other, or (2) has two aggregates 
lying on the opposite side across the cell center such that 
only one independent projection of the chemical gradient 
can be detected. For such a special layout, the Fisher 
information matrix of Eq. (JTTJ) becomes degenerate and 
the cell cannot resolve either gradient steepness p or the 
direction <j>. For M — 2, cr^ diverges at either end, 7Z e = 
or 1, where the system has a single independent sensor 
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FIG. 4. (color online). Sensing limit under varying receptor 
aggregate layout and cell orientation (by non-adaptor recep- 
tors), (a) dp vs. TZe (n = 1) at M = 2 (dashed dot), 3 
(solid) and 4 (dashed). R e was partitioned into 100 equis- 
paced bins within the interval (0, 1) and <Tp was the geomet- 

(23) 

ric mean over 10,000 instances in each bin. Standard errors 
around the means are shown at several R e values, (b) 
(lower 4 curves) and (upper 4 curves), equidistance layout 
(solid) vs. two anisotropic layouts without receptor coopera- 
tivity: (i) Tie = 0.38, 0i = 0.156tt, 2 = 1.08tt and 3 = 1.99tt 
(dashed) and (ii) Tie = 0.49, 0i = 0.513tt, 2 = 1.66tt and 
0s — 1.68-7T (dashed dot), or with receptor coupling of size 
n = 10 for layout (ii) (dotted). Cell rotation is simulated by 
relatively varying cf> from to it. and o\ were calculated 



by Eq. (JSJ) with parameter values iVtot = 80,000, M 
p = 0.1, c = 0.1, a = 0.1, P = 1, and AE = 0. 
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[Fig. Ha)]. For M = 3 or M = 4, cr^ diverges at K e = 1 
and a singularity of the sensing limit may also happen 
at TZ e = 1/2 or 2/3, respectively. Notice that, as the 
number of receptor aggregates M increases, the extent of 
singularity decreases because receptor aggregate layouts 
with the corresponding value of TZ e become less likely to 
be degenerate by chance. 

In general, one can verify that may become sin- 
gular at TZ e = (M — 2)/(M — 1) and lZ e = 1 because 
some receptor aggregates locate so close to each other 
that they effectively sense a same local ligand concentra- 
tion. Suppose among all immediately neighboring aggre- 
gates along the cell circle there exist k angles that are 
equal to or greater than the average angle 2ir/M (i.e., 
the rest M — k angles are smaller than 2ir/M). Let 
be the sum of these k angles. We rewrite the anisotropic 
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factor of Eq. ((23J) as: 



n e = 



9 - 27rfc/Af 
2tt(M - 1)/M 



(24) 



which changes from to 1 as ranges from 2nk/M to 
27r. The aggregate layout can be classified into three 
cases when the value of k changes (see Fig. [5^b) for ex- 
ample layouts when M = 3). (i) When k = 1 and 
6 = 27r, lZ e = 1, the system has a single effective ag- 
gregate at one location (lZ e = 1), and therefore cannot 
resolve parameters p or (f>. The system is non-degenerate 
as long as 8 ^ 2ir. (ii) When k — 2 and 8 = 2ir, 
lZ e = (M — 2)/(M — 1), there exist two effective aggre- 
gates. In the special case in which these two effective ag- 
gregates locate at two opposite sides on the circle across 
the cell center, the Fisher information matrix degener- 
ates. The system is non-degenerate as long as 8 ^ 2tt. 
(iii) When k > 3, no degeneration appears. As an illus- 
tration for M = 3, FiglSfb) shows that a singularity in 
<7p may happen at 1Z e = 0.5 (corresponding to the last 
layout in Fig. EJa)). However, the majority of layouts 
with lZ e = 0.5 are non-degenerate. For instance, all lay- 
outs with k = 1 and G = 4ir/3 has lZ e = 0.5 but none of 
them is degenerate. 

Both ct 2 and c 2 , are sensitive to the cell orientation 
under anisotropic aggregate layout [see Fig. @|b) for a 
non-adaptive case]. The effect of receptor cooperativity 
is preserved under anisotropic receptor aggregate layouts 
because £„ is geometrically independent under a shallow 
gradient. By changing the cell orientation, improvement 
in estimating one parameter of p or </> is gained at the ex- 
pense of the other. This result coincides with the study 
by Hu et al. who showed that an elliptic cell cannot 
simultaneously improve limits c 2 and by elongating 
the cell body. As an optimal configuration, the isotropic 
layout of receptor aggregates improves sensing both p and 
<\> in most cases and is insensitive to cell orientation as 
indicated by Eq. ([12")) and as shown in Fig.@Jb). This re- 
sult recapitulates the one derived by Berg and Purcell |l| 
who showed that uniformly distributed receptor patches 
on the membrane maximize independent ligand influx to 
the cell and thus sensitizes the concentration measure- 
ment. However, as shown by Fig.[4|b), proper cell orien- 
tation can improve estimation of one of the parameters 
(p or (/)), which might be a desirable task if sensing one 
parameter is more critical than the other. For example, 
accurate directional sensing may be more important than 
steepness sensing in chemotaxis so that a cell can iden- 
tify swimming direction more efficiently. For a cell with 
anisotropic layout of receptor aggregates, routine reori- 
entation is required to improve sensing performance (e.g., 
tumbling of E. coli) . 

The above results of minimal number and optimal lay- 
out of receptor aggregates were obtained under the as- 
sumption that the background ligand concentration cq is 
known. One can also consider Cq as an extra parame- 
ter to be estimated by the cell. In this case, merely by 
counting the number of parameters (p, 4> and cq) the cell 



Regular 
layout 



Re=0 



Receptor 
aggregate 




0.2 0.4 0.6 0.8 

Anisotropic factor (R e ) 



FIG. 5. (a) Sample layouts of receptor aggregates (M = 3). 
Note that lZ e — 0.5 corresponds to any scenario of the cell 
with two effective aggregates on cell surface and only the one 
shown with two aggregates on the opposite sides across the 
cell center is degenerate. Dashed lines inside the cell circle 
give a reference to an equiangular partition, (b) Correspon- 
dence between possible k values and lZ e and their relationship 
with <Tp (shown as the geometric mean with s.e.m.). 



needs in minimum three independent receptor aggregates 
as concentration sensors to reliably reconstruct the gra- 
dient. To compute the CRLB's, one must also calculate 
the derivative dc m /dco and obtain a 3-by-3 Fisher infor- 
mation matrix. The primary results about the sensing 
limits remain little changed. It is straightforward to ver- 
ify that under isotropic receptor layout the variances of 
parameter estimates: 



(2 + /3V)(l + c ) 2 (l + aCo )- 
AW^l -a) 2 c 2 



(l+/3V/2)p 2 ' 
(l + c ) 2 (l+ac ) 5 



2 _ 

A t nt£„(l 



ay 



(25) 
(26) 
(27) 



cC approaches the 



Under a shallow gradient ((3p 1), cr 2 
one derived in Eq. (|12p . while the variance of the direc- 
tional sensing is unchanged. The variance <7 2 is in a 
similarly form of Eq. ([7]) for detecting local concentration 
c m by a single receptor aggregate. 
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IV. DISCUSSION AND CONCLUSION 

It is well-known that adaptation brings the receptor 
activity back to pre stimulus level, allowing the system 
to remain sensitive to chemical changes in the future. 
This mechanism enables the system to respond within 
a wide dynamic range of chemical concentration. Here 
we examine the role of adaptation in terms of chemical 
sensing accuracy, especially when presented with small 
signals or small changes in the signal, where stochastic 
noise may become overwhelming. Our study showed that 
receptor cooperativity improves sensing accuracy only 
within a limited background concentration range if the 
receptor aggregates do not undergo adaptation. By con- 
trast, receptor adaptation, especially precise or near pre- 
cise adaptation, maintains the receptor sensory system to 
operate within the parametric region where receptor co- 
operativity is beneficial. As to our knowledge, the finding 
in this work for the first time connects receptor adapta- 
tion and cooperativity to noise filtering in chemical sens- 
ing. 

The sensing limit of Eq. (fT2|) was derived from the 
uncertainty of an instantaneous receptor sampling of a 
gradient field. The sensing accuracy can be improved by 
a cell integrating independent receptor state configura- 
tions over time The time interval between two inde- 
pendent measurements is determined by the correlation 
time r c that accounts for relaxation times of ligand dif- 
fusion, ligand-receptor binding [5| and receptor dynam- 
ics. Because the amount of time t available for averaging 
is limited by dynamics of the intracellular signaling cir- 
cuit, the number of independent measurements is about 
[13 ■ Therefore, the temporal averaging reduces the 
instantaneous limit to ((Tp)t « 2r c ap/t when t 3> 2r c . 
Ligand concentration is encoded as receptor occupancy 
by the ligand-receptor binding and is then transduced 
into receptor activity by the MWC model. Stochastic 
dynamics of each of these signal transduction stages con- 
tributes additional noise and thus increases the corre- 
lation time r c , which consequently elevates the sensing 
uncertainty within the given time frame t. 

Our study does not consider the temporal dynamics 
of ligand-receptor binding and receptor state switching. 
Skoge et al. [Tj recently showed that receptor cooperativ- 
ity may significantly slow down receptor state switching 
(described by Glauber dynamics [2l[ originally developed 
for studying time dependence of the Ising model) to an 
extent such that the system cannot effectively relay the 
information of ligand concentration changes in time [r c 
becomes comparable or greater than the averaging time 
t]. However, there is still a lack of direct experimental 
evidence on how receptor coupling quantitatively mod- 
ulates receptor dynamics. Our results reveal that the 
benefit provided by receptor coupling to chemical sens- 
ing by instantaneous receptor activity configurations may 
compensate the potentially detrimental effect of receptor 
activity slowdown, in particular by adapted receptor ar- 
rays. 



Another possible source of noise at the ligand-receptor 
interaction level is due to ligand rebinding when a disso- 
ciated ligand diffuses back again to bind onto the receptor 
aggregate before it escapes into the bulk medium. The 
extent of ligand rebinding and its effects on the accuracy 
of ligand sensing depends on the size and density of a 
receptor aggregate [HI . Theory [23| and simulation [24j 
showed that ligand rebinding does not change the equi- 
librium occupancy of individual receptors and therefore 
does not affect the sensing limits obtained from instanta- 
neous measurement. However, significant ligand rebind- 
ing does introduce extra fluctuations by decreasing the 
effective rate constants for ligand association with and 
dissociation from the receptor, which increases the corre- 
lation time r c and thus reduces the number of indepen- 
dent measurements within a fixed averaging time. 

Chemoreceptor cooperativity in bacteria was predicted 
by Ising models and MWC models used for analyzing 
dose response data in E. coli. More recently receptor cou- 
pling is evident by the hexagonal organization of trimer- 
of- Aimers of chemoreceptors resolved by electron cryoto- 
mography as a conserved architecture in a wide variety 
of bacteria [2||. The spatial organization of eukaryotic 
chemoreceptors remains to be fully resolved, even though 
higher-order receptor arrays in chemotactic eukaryotes 
have also been observed (26[, suggesting possible recep- 
tor coupling. Our model of the receptor organization on 
the cell surface is general and can be parameterized to 
study systems of coupling or non-coupling receptors with 
or without adaptation. 

Berg and Purcell [l| showed that uniformly distributed 
receptor patches over the cell surface is optimal for re- 
ducing interactions between nearby receptors and there- 
fore maximizes ligand intake by the cell, where ligands 
were considered spatially homogeneous in the environ- 
ment. Here, we consider the detection of spatial asym- 
metry of ligand concentrations around the cell due to the 
chemical gradient. By examining the effect of geomet- 
ric layout of receptor aggregates on the cell surface, we 
showed that anisotropic receptor aggregate distribution 
generates a trade-off between the sensing limit of gra- 
dient steepness p and that of the gradient direction (j> 
(Figj4j) . The isotropic layout represents an optimal con- 
figuration, in which the sensing accuracy is insensitive to 
cell orientation. 

In summary, we study the effect of receptor coopera- 
tivity and adaptation on the chemical sensing limit by 
evaluating the Cramer-Rao lower bounds from the in- 
stantaneous global state of receptor activity. Our results 
showed that receptor cooperativity with receptor adapta- 
tion increases gradient sensing accuracy (by lowering the 
CRLB) for small signal changes across a wide dynamic 
range of background concentration. This result is also ap- 
plicable to concentration sensing at a single aggregate lo- 
cation (Eq. (JT))). It remains largely unknown whether or 
how a chemotactic cell achieves its sensing limit. Recep- 
tor internalization [271 ] might improve chemical sensing 
by helping the cell membrane to function as an absorbing 
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surface, an ideal device that operates at the fundamen- 
tal Berg-Purcell limit. Maximum likelihood estimation 
(MLE) from ligand-receptor bin ding time series [28[ or 
instantaneous receptor states 0, Il9| was also suggested 
as a possible approach to the physical limit. The answer 
to how a cell mechanistically integrates information in 
time and space remains speculative and requires further 



experimental and theoretical investigations. 
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